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Ekpyrotic models and their cyclic extensions solve the standard cosmological flatness, hori- 
zon and homogeneity puzzles by postulating a slowly contracting phase of the universe prior 
to the big bang. This ekpyrotic phase also manages to produce a nearly scale-invariant spec- 
trum of scalar density fluctuations, but, crucially, with significant non-gaussian corrections. 
In fact, some versions of ekpyrosis are on the borderline of being ruled out by observations, 
while, interestingly, the best-motivated models predict levels of non-gaussianity that will be 
measurable by near-future experiments. Here, we review these predictions in detail, and 
comment on their implications. 
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I. MOTIVATION AND INTRODUCTION 

The standard big bang cosmology is hugely successful in describing the evolution of our universe 
from the time of nucleosynthesis onwards. However, a central assumption is that the universe 
started out in a hot big bang, and in a special state: extrapolating back from current knowledge, 
we know that early on the universe must have been very flat, homogeneous and isotropic, with in 
addition small density perturbations with a nearly scale-invariant spectrum and a nearly gaussian 
distribution. Hence, the "initial" state of the universe was far from random, and its specialness 
prompts us to try and explain it via a dynamical mechanism. 

The most studied such mechanism is the model of inflation, which assumes that there was a 
phase of rapid, accelerated expansion preceding the hot big bang; for a comprehensive review see [1]. 
Such a phase can be modeled by having a scalar field (the "inflaton") with a positive and suitably 
flat potential. Inflation has the property of flattening the universe, so that, if it lasts long enough, 
the flatness of the "initial" state can be explained. Moreover, inflation possesses the remarkable 
byproduct that it generates nearly scale-invariant spectra of scalar and tensor perturbations by 
amplifying quantum fluctuations. The predicted scalar perturbations are in excellent agreement 
with current observations, but the tensor perturbations have yet to be observed - their discovery 
would be a strong indication for the correctness of the inflationary picture. However, inflation 
also presents a number of conceptual problems: for example, even though the inflationary phase is 
supposed to erase all memory of initial conditions, this is not really the case. In order for inflation 
to start in a given patch of space, that patch must be reasonably smooth over several Planck 
lengths and the inflaton field must have a small initial velocity (the "patch" and "overshoot" 
problems respectively, see e.g. [1, 2]). Also, it has been realized not long ago that inflation is 
geodesically incomplete towards the past, which means that the predictions of the theory depend 
on the specification of data on a spacelike initial hypersurface [3]. In other words, inflation requires 
its own initial conditions. Hence, if inflation is correct, it will only form a part of the story. More 
worrying is the problem of unpredictability, which is associated with the quantum nature and the 
effectiveness of inflation. Inflation ends when the inflaton field oscillates around a minimum of 
its potential, and "reheats" the universe by decaying into standard model particles. However, for 
generic initial conditions there will always be regions in which rare but large quantum fluctuations 
kick the inflaton field back up its potential and keep a fraction of the universe in the inflationary 
phase. In most of the concrete realizations of inflation, the region that keeps inflating expands 
so fast that it quickly dominates the overall volume of the universe. Hence, inflation never ends 



3 



and the global picture of this process of "eternal inflation" is that of an empty de Sitter universe 
punctured by an infinite number of small pockets where inflation has ended (at a random time) 
[4]. Because inflation ends at a random moment in these pocket universes, the pockets might 
have become sufficiently flattened or not, they might have acquired scale-invariant perturbations 
or not. Without a measure which would determine the relative likelihood of the various pockets, 
it becomes difficult to know exactly what eternal inflation predicts! These problems do not mean 
that the idea of inflation is wrong, but, if inflation continues to be supported by observations, they 
will have to be addressed. In the meantime, the seriousness of these open problems means that it 
is worthwhile considering alternative models for the early universe in parallel. 

The present review deals with one such model in particular, namely the ekpyrotic model and its 
extension, the cyclic universe; for a comprehensive overview see [5]. In this model, the inflationary 
phase is replaced by the ekpyrotic phase, which is a slowly contracting phase preceding the big 
bang. The ekpyrotic phase can be modeled by having a scalar field with a negative and steep 
potential. As described in detail below, it also manages to flatten a given region of the universe, 
and generates nearly scale-invariant scalar perturbations, but no observable tensor fluctuations. 
At the linear level, the scalar fluctuations are virtually indistinguishable from the perturbations 
produced by inflation, but at higher orders the predictions differ. Since primordial gravitational 
waves might turn out to be rather elusive to measure over the coming years, the most promising 
way of distinguishing between alternative models of the early universe is therefore by studying 
these higher-order, non-gaussian signatures. 

There is a simply, intuitive argument for why the predictions regarding higher-order corrections 
to the linear perturbations should differ for models of inflation and ekpyrosis. For a scalar field 
fluctuation Sip, the semi-classical probability density is roughly given by e~ SE ^ s<p \ where Se{6<p) is 
the euclidean action [6]. Since inflation requires a very flat potential, the inflaton is an almost free 
field. For a free field, the action is quadratic in the field, and hence the probability distribution 
is simply a gaussian distribution. For an exact gaussian distribution the 3-point function (S(p 3 ) 
vanishes, and hence for inflation, where the field is almost free, we would expect the 3-point function 
to be non-zero, but small. For ekpyrosis, on the other hand, the potential is steep, and hence the 
scalar field is necessarily significantly self-coupled. This has the consequence that ekpyrotic models 
generally predict significant levels of non-gaussianity. In fact, some versions of ekpyrosis are already 
on the borderline of being ruled out by observations, while the best-motivated models predict values 
that are measurable by near-future experiments. Thus, the non-gaussian predictions are crucial 
in assessing the viability of various cosmological models, and promise to significantly enhance our 
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understanding of the physics of the early universe. 

The plan of this review is to start with a brief summary of the main ideas behind ekpyrotic 
and cyclic models of the universe. We will then discuss in some detail the generation of linear 
cosmological perturbations (a good understanding of the linear perturbations greatly facilitates an 
understanding of the higher-order ones), before turning to the main subject of the review, namely 
the non-gaussian corrections to these linear perturbations. We will conclude with a discussion of the 
non-gaussian predictions and in particular their observability and relation to current observational 
limits, as well as the consequences of a potential detection. 



II. EKPYROTIC AND CYCLIC COSMOLOGY 



The ekpyrotic phase is the cornerstone of ekpyrotic and cyclic models of the universe: it is a 
conjectured, slowly contracting phase preceding the big bang, and it resolves the standard cos- 
mological puzzles [7, 8]. The main feature of ekpyrosis is that during this phase the equation of 
state 

w = -»l (1) 
P 

is very large (here p and p denote the average pressure and energy density of the universe). Let 
us briefly explore the most direct consequences of such an ultra-stiff equation of state. Consider a 
Friedmann-Robertson- Walker (FRW) metric 1 

ds i = _ d e + a {tf (y^2 + r2 ^) ( 2 ) 

where a(t) denotes the scale factor of the universe and k = —1,0,1 for an open, flat or closed 
universe respectively. If the universe is filled with a number of fluids interacting only via gravity 
and with energy densities pi and constant equations of state Wi , then the equations of continuity 

Pi + 3-( Pi +Pi) = (3) 
a 

(where dots denote derivatives with respect to time t) imply that they will evolve according to 

Pi oca~ 3{1+w >\ (4) 

The Einstein equations for this system contain a constraint equation, better known as the Fried- 
mann equation, which involves the Hubble parameter H = a/a: 

H 2 = U=^ + ^ + ^f + ^f + ... + ^-). (5) 

3 I a 2 a 3 a 4 a e a 3(i+^) / v ; 



1 I will mostly use natural units h = c = 1 and 8tyG = M F 



FIG. 1: The potential during ekpyrosis is negative and steeply falling; it can be modeled by the exponential 
form V{<j>) = -V e- cr t>. 

The pifiS are constants giving the energy densities at scale factor a = 1 of the various constituents 
of the universe: we consider the universe to be composed of non-relativistic matter (subscript m), 
radiation (r) and the energy density associated with anisotropies in the curvature of the universe 
(a). In addition, we consider there to be ekpyrotic (scalar) matter, denoted by the subscript 0, 
and, as usual, there is a contribution due to the average curvature of space. 

As the universe contracts, components whose energy density scales with a higher negative 
power of the scale factor a will successively come to dominate, first matter, then radiation, then 
anisotropies and eventually, since ^ 1 by assumption, the ekpyrotic matter. This means 
that the relative energy densities in curvature and anisotropies, for example, become smaller and 
smaller, the longer the ekpyrotic contracting phase lasts. In other words, if ekpyrosis lasts long 
enough, the flatness problem is solved. We will make this statement quantitative below. Strictly 
speaking, for the flatness problem to be solved, all we need is a matter component with w > 1. In 
the next section, we will see that for realistic ekpyrotic models, typically 3> 1. In passing, we 
should also point out that there is no horizon problem in ekpyrotic and cyclic models, as there is 
plenty of time before the big bang for different parts of our currently observable universe to have 
been in causal contact with each other. 

But what form of matter can have the large equation of state that we require? A simple way 
to model the ekpyrotic matter is to have a scalar field <j> with a steep and negative potential V {(/>)■ 
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A concrete example is provided by the negative exponential 

V{fa) = -V e~ c t (6) 

where Vq and c are constants - see Fig. 1. In the context of string theory, such scalar fields appear 
very naturally, and the ekpyrotic potential can then correspond to an attractive force between 
branes - this picture will be briefly described below. 

Given an explicit form of the potential, such as (6), we can solve for the evolution of the universe. 
In fact it is straightforward to generalize the treatment to having many scalars fa with potentials 
Vi(fa). Then, in a flat FRW background and neglecting other matter components, the equations 
of motion become 

fa + 3Hfa + Vi^ = (7) 

and 



H 



2 



(8) 



where = (dVi/dfa) with no summation implied. If all the fields have negative exponential 
potentials Vi(fa) = — V^e _c< ^ and if Cj S> 1 for all i, then the Einstein-scalar equations admit the 
scaling solution 

a = {-tfl\ fa = l H -^J~2t), ~ = Y.\- ( 9 ) 
Thus, we have a very slowly contracting universe with (constant) equation of state 

, = aft'-w.) * t>1 . (10) 

We are using a coordinate system in which the big crunch occurs at t = 0; in other words, the 
time coordinate is negative during the ekpyrotic phase. Here, the parameter e corresponds to the 
fast-roll parameter and is typically of 0(100); its definition is identical with that in inflation, where 
its value is typically of 0{ 1/100) and where, correspondingly, it is called the slow-roll parameter. 

Using this explicit solution, we can get an idea for how long the ekpyrotic phase has to last 
in order for the flatness problem to be solved. Quantitatively, the problem can be formulated as 
follows: dividing the Friedmann equation (5) by H 2 we can see that the fractional energy density 
stored in the average curvature of the universe is given by 
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At the present time, observations imply that this quantity is smaller than 1CT 2 in magnitude [9]. 
If we assume a radiation-dominated universe, which is a good approximation for this calculation, 
then aH oc t -1 / 2 and hence, if we extrapolate back to the Planck time, the fractional energy density 
in curvature must have been smaller than 

^10- 2 « 10" 62 , (12) 
to 

an incredibly small number. However, from (9), we can see that during the ekpyrotic phase the 
scale factor a remains almost constant, while the Hubble parameter H oc t -1 . Hence aH grows by 
a factor of 10 30 as long as 

\tek— beg\ ^ e l^efc— end\t (13) 

where the subscripts ek — beg and ek — end refer to the beginning and the end of the ekpyrotic 
phase respectively. As will be discussed in the next section, we need t e k- e nd ~ — 10 3 Mpl in order 
to obtain the observed amplitude of cosmological perturbations, so that we need 

\t ek . beg \ > W 33 M p l = W- 10 s. (14) 

This is the minimum time the ekpyrotic phase has to last in order to solve the flatness problem. 
Cosmologically speaking, this is a very short time, attesting to the effectiveness of the ekpyrotic 
phase. 

Before discussing the cosmological perturbations produced during the ekpyrotic phase, it is 
useful to provide a quick overview of how the ekpyrotic phase might fit into a more complete 
cosmological model. The crucial ingredient in any such model is the proposed mechanism for how 
the ekpyrotic contracting phase (with H < 0) and the subsequent radiation-dominated expanding 
phases (with H > 0) should link up. The Einstein equations provide the relation 

H = -\{p + p). (15) 

All forms of matter that are currently known to exist obey the null energy condition 

p + p^0 (NEC), (16) 

which implies H ^ and which thus precludes a smooth transition between a contracting and an 
expanding universe. This leaves two possibilities for achieving such a transition: either the NEC 
is violated during the transition, or the transition is classically singular. 

In new ekpyrotic models [10-12], a smooth reversal from contraction to expansion is achieved 
by adding a further matter component to the universe which can violate the NEC. The particular 
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example that these models consider is the so-called ghost condensate, which corresponds to the 
gravitational equivalent of a Higgs phase [13]. It is not clear yet whether or not the ghost condensate 
can be obtained from a fundamental theory such as string theory [14] (in the more restricted 
framework of quantum field theory it seems impossible to construct a stable ghost condensate model 
[15]); however, it is interesting that string theory contains many objects (orientifolds, negative- 
tension branes) which do violate the NEC. Of course, simply adding such a component is not 
enough: it must become relevant as the universe contracts, and vanish again as the universe 
expands. The simplest way in which to achieve this is by assuming that the ghost condensate itself 
also plays the role of the ekpyrotic matter, and that after the transition to expansion, it decays into 
ordinary matter fields. This scenario requires the ghost condensate to possess both a particular 
form for its kinetic term and a particular potential; for details regarding possible realizations see 
[10, 12]. 

The cyclic model of the universe [16, 17] is based on the braneworld picture of the universe, 
in which spacetime is effectively 5-dimensional, but with one dimension not extending indefinitely, 
but being a line segment, see Fig. 2. The endpoints of this line segment (orbifold) are two (3 + 1)- 
dimensional boundary branes. In the full string theory setup, there is in addition a 6-dimensional 
internal manifold at each point in the 5-dimensional spacetime, for a total of 11 dimensions [18]. 
This description of the universe stems from string theory, and in particular the duality, known as 
Hofava-Witten theory [19], between 11-dimensional supergravity and the Eg x Eg heterotic string 
theory. All matter and forces, except for gravity, are localized on the branes, while gravity can 
propagate in the whole spacetime. Our universe, as we see it, is identified with one of the boundary 
branes and, as long as the branes are far apart, can interact with the other brane only via gravity. 
The cyclic model assumes that there is an attractive force between the two branes, which causes the 
branes to approach each other. This force is modeled by a potential of the form shown in Fig. 3. 
Note that the potential incorporates an ekpyrotic part. From the higher-dimensional point of view, 
the ekpyrotic phase has the rather non-intuitive property that it flattens the branes to a very high 
degree. Eventually the two branes collide and bounce off each other. It is this collision that, from 
the point of view of someone living on one of the branes, looks like the big bang. Classically, the 
collision is singular, as the orbifold dimension shrinks to zero size. The collision is slightly inelastic 
and produces matter and radiation on the branes, where the standard cosmological evolution now 
takes place. However, due to quantum fluctuations, the branes are slightly rippled and do not 
collide everywhere at exactly the same time. In some places, the branes collide slightly earlier, 
which means that the universe has a little bit more time to expand and cool. In other places, 
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orbifold dimension 



FIG. 2: The braneworld picture of our universe. Think of a sandwich: the 5-dimensinonal bulk spacetime is 
bounded by two 4-dimensional boundary branes. There is no space "outside" of the sandwich, but the branes 
can be infinite in all directions perpendicular to the line segment (orbifold). In the M-theory embedding, 
there arc 6 additional internal dimensions at each point of the sandwich. 

the collision takes place slightly later, and those regions remain a little hotter. This provides a 
heuristic picture of the way temperature fluctuations are naturally produced within the model. 
Shortly after the branes have separated, the distance between the boundary branes gets almost 
stabilized, but the branes start attracting each other again very slightly. This very slight attraction 
acts as quintessence, and is identified with the dark energy observed in the universe. After a long 
time, and as the branes become closer again, they start attracting each other more strongly so that 
we get another ekpyrotic phase and eventually another brane collision with the creation of new 
matter. In this way, a cyclic model of the universe emerges. Before continuing, we should mention 
the main open issues related to the cyclic model: the first one concerns the potential, which at 
this point is simply conjectured. It will be important to see if a potential of the required shape 
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FIG. 3: The potential for the cyclic universe integrates the ekpyrotic part and a quintessence epoch, but 
is irrelevant at the brane collision. A possible form for the potential is V{<j>) — Vb(e h ^ — e~ c ^)F((f)), with 
t«l,c>l and F((f>) tends to unity for <f> > <j) en d and to zero for (f> < <j) en d- Reproduced with permission 
from [8], 

can be derived from micro-physics. And the second is the brane collision, which so far has been 
extensively studied at the classical and semi-classical level [20], but a full quantum treatment has 
remained elusive. 

In the discussion above, we have mostly focussed on models involving one effective scalar field. 
However, there are two good reasons to extend the analysis to two or more scalars: first, in 
embedding the ekpyrotic and cyclic models in M-theory, there are two universal scalars, namely 
the radion mode (which determines the distance between the branes) and the volume modulus 
of the internal 6-dimensional manifold [21]. There can be many more scalar fields (such as the 
shape moduli of the internal space), but we always must consider these two universal scalars. And 
secondly, as we will see in the next section, it is much more natural to generate a nearly scale- 
invariant spectrum of curvature perturbations (in agreement with observations) in models with 
two scalars than in models with only one. However, multi-field ekpyrotic models present some 
qualitatively new features, which we discuss briefly here. 
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The 4-dimensional effective action 
S = 



1 



g[R - ^(dfaY - ^{d<t> 2 f - V(0i, 2 )] 



(17) 



can be obtained as the low-energy limit of Hof ava-Witten theory, where 0i and 02 are related by a 
field redefinition to the radion and the internal volume modulus [22]. We are assuming that during 
the ekpyrotic phase, both fields feel an ekpyrotic-type potential, e.g. 



V(<h,<h) 



-Vie- 01 * 1 - V 2 e~ 02 ^ 2 . 



(18) 



Then it is much more natural to discuss the dynamics in terms of the new variables a and s 
pointing transverse and perpendicular to the field velocity respectively [23, 24]; they are defined, 
up to unimportant additive constants which we will fix below, via 

h<j>i + 0202 



a 



a 



>W2 - 9291 
& 



(19) 



with & = (<pf + 02 ) 1 ^ 2 - It is a l so useful to define the angle 9 of the trajectory in field space, via [25] 

(20) 



cos I 



sm( 



a 



a 




FIG. 4: After a rotation in field space, the two-field ekpyrotic potential can be viewed as composed of an 
ekpyrotic direction (a) and a transverse tachyonic direction (s). The ekpyrotic scaling solution corresponds 
to motion along the ridge of the potential. Perturbations along the direction of the trajectory are adi- 
abatic/curvature perturbations, while perturbations transverse to the trajectory are entropy/isocurvature 
perturbations. 
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In terms of these new variables, the potential can be re-expressed as 

V ek = -V e^[l + es* + | e 3 /V + ^eV + •••], (21) 

where for exact exponentials of the form (18), one has K3 = 2\^2(cf — c\)j\c\C2\ and K4 = 4(c^ + 
c i) / ( c i c l ( c i + c i))- However, in the absence of a microphysical derivation of the potential, we 
will simply take ^3,^4 ~ 0(1) and express all results in terms of ^3,^4. See also Fig. 4 for an 
illustration of the potential. The ekpyrotic scaling solution becomes 

a(t) = (-t) 1/t a = -^\n(-^eVot^j s = 0, (22) 

with the angle 6 being constant. The solution corresponds to motion along a ridge in the potential, 
as is evident from the figure. Hence, in contrast to the single field case, the multi-field ekpyrotic 
background evolution is unstable to small perturbations [26, 27]. This implies that the trajectory 
must be localized near the ridge with extreme precision at the beginning of the ekpyrotic phase, 
the condition being that the field should stray no more than a value of e" 60 (at best) in Planck 
units from the ridge at the beginning of ekpyrosis [11]. Thus, at first sight, it looks as if the multi- 
field ekpyrotic phase has not managed to solve the problem of initial conditions. However, there 
currently exist two approaches addressing this issue: the authors of [11] considered the existence of a 
"pre-ekpyrotic" phase during which the potential is curved upwards and during which the trajectory 
is localized. Meanwhile, in the context of the cyclic universe, there is a natural resolution of the 
issue of initial conditions, not involving any new ingredients of the model: indeed, the multi-field 
cyclic universe selects those regions that happen to correspond to trajectories sufficiently close to 
the ridge, in the sense that these regions are vastly amplified over the course of one cycle due to 
the phases of radiation, matter and dark energy domination (note that the ekpyrotic phase shrinks 
the universe by a negligible amount). At the same time, the regions corresponding to trajectories 
not sufficiently close to the ridge (this would include the vast majority of trajectories) do not 
undergo a full ekpyrotic phase, and after these regions undergo chaotic mixmaster behavior close 
to the big crunch, they simply collapse (presumably they will end up forming black holes) and stop 
both growing and cycling. In this way the global structure of the universe becomes of the phoenix 
type, in which vast habitable regions are interspersed with small collapsed ones. The important 
point is that the habitable regions, which are the only regions of interest to us here, automatically 
correspond to the regions that had the right "initial conditions" at the beginning of their preceding 
ekpyrotic phase. This is discussed in detail in [28]; see also the essay [29]. 
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III. LINEAR PERTURBATIONS 
A. Single Scalar Field 

In the last section, we have dealt with the classical evolution during the ekpyrotic phase. We will 
now add quantum fluctuations and we will see that, just as in inflation, the quantum fluctuations get 
amplified into classical density perturbations. Hence, on top of resolving the standard cosmological 
puzzles, the ekpyrotic phase can also be the source of the primordial temperature fluctuations 
whose imprint is seen in maps of the cosmic microwave background, provided that the amplitude 
and spectrum of the fluctuations match observations. 

From the study of inflationary models, we have developed the intuition that quantum fluctua- 
tions that get stretched to super-horizon scales turn into classical perturbations, roughly speaking 
because the fluctuations go out of causal contact with themselves, do not remember locally that 
they are in fact fluctuations and end up as local, classical perturbations to the background evolu- 
tion. In inflation, this effect occurs because the horizon is approximately constant in size while the 
wavelengths of the quantum modes get stretched exponentially with time (the scale factor of the 
universe grows exponentially). For ekpyrosis, the scaling solution (9) shows that the scale factor is 
almost constant, so that the mode wavelengths remain almost constant too. However, the horizon, 
which is proportional to l/H ~ t, shrinks rapidly as t — > and hence the modes automatically 
become of super-horizon size 2 . We will now discuss in some detail what amplitude and spectrum 
these modes obtain. We will first concentrate on the single field case, before discussing two fields. 

Since the scale factor evolves very little during the ekpyrotic phase, one is tempted to simply 
turn gravity off as a first approximation, and to consider the theory consisting only of a scalar field 
with a steep and negative potential [32]: 

S = J rf 4 x[-^(%) 2 + V ie ^H (23) 

Then, if we define scalar fluctuations 5(p via <f>i = 4>i(t) + 5(p(t,x), where (pi = ln(— y^V\/2t) 
denotes the background evolution, the equation of motion for the fluctuations is given by 

8<f>- V 2 5<p + V Ml 54> = $, (24) 

2 Since tensor modes/gravitational waves depend on the evolution of the scale factor alone, and since the scale 
factor shrinks imperceptively slowly during ekpyrosis, there are no substantial gravity waves produced during the 
ekpyrotic phase (the background spacetime is almost Minkowski!) [30]. In fact, the dominant gravitational waves 
that are produced from ekpyrosis are those that arise from the backreaction of the scalar fluctuations onto the 
metric, at second order in perturbation theory [31]. 
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where V l p 1( j >1 = —2/t 2 . We then expand the fluctuation field 5(f) into Fourier modes 

d k ^a kX ke ik '* + h.c. (25) 



(2tt) 2 

where the XfcS are the positive frequency mode functions (due to the assumed cosmological symme- 
tries, they depend only on the magnitude k = |k|). We proceed to quantize the field by imposing 
the canonical commutation relations 

[ok, aw] = [at, al] = 0, [a k , a*,] = (2nf5(k - k') (26) 

In the process, the a^s have been promoted to (annihilation) operators, and the vacuum state |0) 
is defined by dk|0) = 0. The mode functions obey the equation of motion 

Xk + k 2 Xk ~ ^Xk = 0, (27) 
which admits the two solutions \k °c e~ lkt {l — -j^),e lkt (l + ^). However, as t — > -co the modes 



should asymptote to the Minkowski space free particle state Xk — > e lkt /V2k (note that in that 
limit (27) reduces to the equation of a simple harmonic oscillator), and this fixes the solution to 
be 

Towards the end of the ekpyrotic phase, we have \kt\ <C 1, and then the solution can be well 
approximated by 

Xk * vSm- (29) 

The quantum fluctuations have a mean that is zero, (0|<5</>|0) = 0. However, the variance A^(/c), 
which is defined by (O|(50 2 |O) = J ^A 2 (fc), does not vanish. It is conventional to write the variance 
as 

Al(k) = Al(k )(^-T°- 1 , (30) 

where fco denotes a reference scale and n s is the spectral index. 

A related concept in momentum space is the power spectrum P(k), defined by 

9-7T 2 

P(k) = \xk\ 2 = ^A 2 (fe). (31) 
It is the Fourier transform of the 2-point correlation function, and we can equivalently define it as 

(CkCk') = fP(k)5 3 (k + k') (32) 
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where isotropy dictates that P only depends on k = |k|. We will find this definition useful later on. 
In our case, we have that at late times 

r f d 3 k f Ank 2 dk 1 , , 

'w^lw^'/ww 1 (33) 

so that the variance is given by 

Aj(*) = ^ n s = l. (34) 

The variance is independent of k, and hence we obtain a scale-invariant spectrum for 5(f). This 
looks very promising! However, we really must include gravity in our analysis, and calculate the 
spectrum for the curvature perturbation £, which is the quantity that is measured to have a nearly 
scale-invariant spectrum of perturbations. 

Once we add gravity, it is easiest to perform the calculation in so-called £-gauge, where the 
perturbations in the scalar field are gauged away and all perturbations are expressed via dilatations 
of the 3-metric: 

5(f) = (35) 
ds 2 = -dt 2 + a 2 (t)e 2 ^ x) dx j dx j , (36) 

where j = 1,2,3. Then, using the background scaling solution (9), the action reduces to an action 
for £ which is given by [33] 



/ 



eg^CduC (37) 



During ekpyrosis, e is typically nearly constant. In fact, in the scaling solution used above, we have 
already made the approximation that e is constant, and with this approximation, the equation of 
motion for £ resulting from the action above is particulary simple: in Fourier space it is given by 

C* + SHCk + ^Ck = 0. (38) 

If we use conformal time r, defined via dt = adr, and the notation ' = jj:, the above equation 
becomes 

+ 2^C' k + k\ k = 0. (39) 

After a further change of variables to y = aC/\/—kr and x = —kr, the equation turns into a Bessel 
equation x 2 ^ + + ( x2 ~ a<2 )y = 0) with a = y 7 r 2 a" /a + 1/4 « 1/2 since a « constant. Hence 
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the solutions are given by the Hankel functions y oc Hy' 2 {— kr), Hy' 2 (— kr) and with the boundary 
condition that we want £ — > e~ %kT j\p2k as r — > -co, we obtain the solution (up to a phase) 3 

C = ^H^i-kr). (40) 

At late times, the variance becomes 

and hence we get a spectral index n s = 3. This spectrum is blue, as there is more power on 
smaller scales, and it is in disagreement with observations [32-34]. Hence, the scale-invariant 
spectrum of the scalar perturbation in the no-gravity theory did not get transferred to the curvature 
perturbation A closer analysis reveals that these two perturbations correspond to two physically 
distinct modes, the former being a time-delay mode to the big crunch, and the latter a local 
dilatation in space. In a contracting universe, these two modes are distinct, and they do not mix. 
It is conceivable that they might mix at the big crunch/big bang transition [35, 36], in which case 
the scale-invariant contribution would be the dominant one on the large scales of interest, but this 
possibility is still insufficiently understood to make definite predictions. As we will show next, this 
is also unnecessary, as there is a very natural entropic mechanism which generates scale-invariant 
curvature perturbations before the big bang, as long as there is more than one scalar field present 4 . 
But before continuing, it might be useful to add a few remarks concerning the validity of our 
approach: indeed, the reader might be worried about the validity of perturbation theory, since the 
background quantities, such as the Hubble rate, as well as the perturbations themselves blow up 
as t — > 0. However, as shown in [33], after switching to synchronous gauge, it is straightforward to 
see that the universe evolves to become closer and closer to the unperturbed background solution, 
and hence perturbation theory is valid. Also, even though the background quantities blow up as 
seen from the 4-dimensional viewpoint, in fact in the higher-dimensional colliding branes picture 
the ekpyrotic phase has the effect of flattening the branes and hence of rendering the curvatures 
small as the big crunch is approached [20]. 

3 Useful asymptotic expressions are H^(x) -¥ yT e i ( I -<"/ 2 -"/ 4 ) when x > a and h£\x) -> -±r(a)(|) Q when 
x <Si a and for a > 0. 

4 Recently, Khoury and Steinhardt have also pointed out that right at the onset of the single-field ekpyrotic phase, 
a range of scale- invariant modes can be produced [37]. However, contrary to the cases that we have discussed so 
far, this adiabatic mechanism requires the universe to already be contracting when the equation of state is still 
near w « — 1. If viable, this mechanism would produce an interesting non-gaussian signal; but as it is currently 
not known how to incorporate this mechanism into a more complete cosmological model, we will not discuss this 
mechanism here. See also [38] for the challenges that this scenario must address. 
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B. Two Fields: the Entropic Mechanism 



As discussed at the end of the last section, it is rather unnatural to consider only a single 
scalar field in the effective theory, since there are two universal scalars that are always present 
in a higher-dimensional context: the radion field, determining the distance between the two end- 
of-the-world branes, and the volume modulus of the internal manifold. But as soon as there is 
more than one scalar field present, one can have entropy, or isocurvature, perturbations, which 
are growing mode perturbations in a collapsing universe [39]. Entropy perturbations can source 
the curvature perturbation, and hence (provided the entropy perturbations acquire a nearly scale- 
invariant spectrum), nearly scale-invariant curvature perturbations can be generated just before 
the bounce [26, 40]. These then turn into growing mode perturbations in the ensuing expanding 
phase. 

For the two scalar fields, we will again assume the potential (21), but with the slight general- 
ization that we allow the fast-roll parameter e to be slowly varying. There are two gauge-invariant 
scalar perturbation modes: the entropy perturbation 5s = cos65(p2 — s'm95(f>i corresponds to per- 
turbations transverse to the background trajectory, see Fig. 4, while the adiabatic, or curvature, 
perturbation ( is the gauge-invariant quantity expressing perturbations along the background tra- 
jectory, see [25, 41] for a detailed exposition. For a straight trajectory (6 = 0), the linearized 
equation of motion for 5s is 



where V ss denotes the second derivative of the potential w.r.t. s. In conformal time, and for the 
re-scaled variable 5S = a(r)5s, we obtain 



To proceed, we must relate a" /a and V ss to the fast-roll parameter e and its derivative w.r.t. 
the number of e-folds of expansion N, where dN = <ilna. By requiring e to vary slowly, what is 
meant is that we will keep terms in de/dN, but not higher-order terms such as d 2 e/dN 2 . Then, by 
differentiating e = a 2 /(2H 2 ) twice, and using a + 3Ha + V a = as well as V ss = V aa , one can 
derive the following expressions, valid to sub-leading order in e: 




(42) 




(43) 




(44) 



(45) 
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Using in addition that aH = (1 + 1/e + e^/e )/(er), Eq. (43) finally reads 

6S " + (k 2 - { 2e J ) 5S = °- (46) 

In analogy with our discussion of the single-field case, this equation can be solved in terms of the 
Hankel functions, supplemented by the boundary condition of approaching the Minkowski vacuum 
state in the far past, to yield (up to a phase) 

SS-^H-V, ,-§d-£ + g>. (47) 
At late times (— kr) — > and we obtain 

5S « -=-i , (48) 

V^(-r)fc" 

implying that at the end of the ekpyrotic phase, the entropy perturbation is given by 

Ss(t ek _ end ) « leVek ^ n k f 2 ■ (49) 

Following the same steps as in the single field case above, it is straightforward to see that the 
spectral index of the entropy perturbation is now given by [26] 

n s - 1 = - - ^. (50) 

The first term on the right-hand side is a gravitational contribution, which, being positive, tends 
to make the spectrum blue. The second term is a non-gravitational contribution, which tends 
to make the spectrum red. A simple way to estimate the natural range of n s is to rewrite the 
above expression in terms of A/", the number of e-folds before the end of the ekpyrotic phase (where 
dM = dln(aH)): 

n - 1 - - - ^ (51) 

In this expression, e(J\f) measures the equation of state during the ekpyrotic phase, which decreases 
from a value much greater than unity to a value of order unity in the last e-folds. If we estimate 
e « N a [42] , then the spectral tilt is 

2 a , 
n s - 1 » — r- - — -.. (52) 

Here we see that the sign of the tilt is sensitive to a. For nearly exponential potentials (a « 1), 
the spectral tilt is n s « 1 + 1/ftf rj 1.02, slightly blue, because the first term dominates. However, 
in the cyclic model the steepness of the potential must decrease in order for the ekpyrotic phase to 



19 



come to an end, and a parameterizes these cases. If a > 1.14, the spectral tilt is red. For example, 
n s = 0.97 for a ps 2. These examples represent the range that can be achieved by the entropic 
mechanism, roughly 

0.97 < n s < 1.02. (53) 

These are in good agreement with the present observational bounds obtained by the WMAP 
satellite, which are n s = 0.96 ± 0.03 at the 2a level [9]. 

Now that we have shown how an approximately scale-invariant spectrum of entropy pertur- 
bations may be generated by scalar fields in a contracting universe, we will discuss how these 
perturbations may be converted to curvature perturbations. Since the entropy perturbations of 
interest are all of super-horizon scales, we can now restrict our study to large scales only, where 
spatial gradients can be neglected. On these scales, the evolution equation for the curvature per- 
turbation is given by [25] 

2 ff /~2~ 

C = — -r-OSs = J -95s. (54) 

Hence, as soon as the background trajectory bends (9 / 0), the entropy perturbations become a 
source for the curvature perturbations. 

There are at least two ways in which such a bending can occur: the first makes use of the 
instability of the two-field ekpyrotic potential, cf. again Fig. 4. If the background trajectory 
strays sufficiently far from the ridge of the potential, the trajectory will turn and fall off one 
of the steep sides of the potential [10, 23, 24]. The turning of the trajectory then immediately 
results in the conversion of entropy into curvature perturbations. Since this conversion occurs 
during the ekpyrotic phase, we will term this process ekpyrotic conversion. It is straightforward 
to estimate the amplitude of the resulting curvature perturbation (its spectrum will be identical 
to the spectrum of the entropy perturbations, as Eq. (54) is /c-independent): if we approximate 
the entropy perturbation as remaining constant during the conversion process, and assume a total 
bending angle of order unity, f 9 ~ 0(1), then the resulting curvature perturbation after conversion 
will be given by 

I 2 

Cconv—end ~ \ $Sek—end- (55) 

V e e fc 

We should mention straight away that the approximations just made will not be good enough in 
calculating the non-gaussian corrections to the linear calculation, but for the present purposes, 
they will do. Since the fast-roll parameter e e u ~ O(10 2 ), we find that 

Cconv—end ~ "Jg ^ s ek~end- (56) 
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FIG. 5: After the ekpyrotic phase, the trajectory in scalar field space enters the kinetic phase and bends - 
this bending is described by the existence of an effective repulsive potential (the potentials are indicated by 
their contour lines). A trajectory adjacent to the background evolution can be characterized by the entropy 
perturbation 6s(t e k- e nd) at the end of the ekpyrotic phase, leading to a corresponding off-set 8s(tb en d), or 
equivalently 5V(tb en d), at the time of bending. 

The second way in which a bending of the trajectory can occur is suggested by the embedding 
of the cyclic model in M-theory [26], and applies when the background field trajectory remains 
straight throughout the ekpyrotic phase. At the end of the ekpyrotic phase, the potential turns off, 
and the evolution becomes dominated by the kinetic energy of the two scalar fields. This kinetic 
phase corresponds to the final approach of the branes in the higher-dimensional picture. During 
this approach, there is a generic effect that occurs, but that cannot be seen in the 4-dimensional 
effective theory. The branes that are bound to collide with each other are of opposite tension. Now, 
it turns out that just before the collision, it always happens that at the location of the negative 
tension brane, the internal 6-dimensional manifold tries to shrink to zero size [21]. However, just 
about any type of matter present on the negative tension brane will smoothly cause the internal 
manifold to grow again [43]. This effect is due to the special properties of gravity on a negative 
tension object. When this effect is translated back into the effective theory that we have been 
using all along, the location in field space where the internal manifold reaches zero size is given by 
the (f>2 = line. This line thus constitutes a boundary to field space. And the presence of matter 
on the negative tension brane gives rise to en effective repulsive potential in the vicinity of the 
4>2 = line. Hence, during the kinetic phase, the background trajectory automatically bends, just 
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before the trajectory shoots off to — oo where the brane collison/big bang occurs. What we have 
just discussed is a concrete example originating from string theory of how a bend in the trajectory 
can occur during the kinetic phase. However, more generally our results will apply whatever the 
microphysical cause of the bending and of the effective repulsive potential, see Fig. 5 for an 
illustration of the general case that we have in mind. Again, it is quite straightforward to estimate 
the amplitude of the curvature perturbation resulting from this process of kinetic conversion: 
On large scales, the linearized equation of motion for the entropy perturbation is given by 

5s + 3H5s + (v ss + 3# 2 ) 5s = 0, (57) 

where, incidentally, we have the useful relation 9 = —V s /&. Then, during the kinetic phase and 
away from the repulsive potential, the Einstein equations immediately yield 

1 -V2 

H = — , a = -= . (kinetic phase) (58) 

ot yj2>t 

Thus, (57) simplifies to 5s + 5s/t = which implies that the entropy perturbation grows logarithmi- 
cally during the potential-free kinetic phase. We can ignore this insignificant growth. However, the 
evolution of the entropy perturbation during the process of conversion turns out to be important. 
We can estimate it by assuming that the trajectory bends with a constant 9 ~ I /At, where At 
denotes the duration of the conversion process. We can further assume that the repulsive potential 
depends only on (p2- Then 9, V s , V ss can all be related to V^ 2 , Vfofo, evaluated during the conver- 
sion, and it is not difficult to show that this leads to V ss ~ <p±9 / {tbendfa) , where tbend corresponds 
to the time halfway through the bending of the trajectory [44]. For the particular example where 
the cyclic model is embedded in M-theory, we have that <\>\ = — V3(t>2, and specializing to this 
example, we have V ss « (2 — 3)9 . Hence, (57) becomes (where we can neglect the term in 5s) 

5s + 69 2 5s « 0, (59) 

and thus, during the conversion, the entropy perturbation evolves sinusoidally 

5s rj cos[uj(t - t conv _ beg )]5s(t ek _ end ) , (60) 

where t con i)-beg denotes the time at which the trajectory starts to bend, and oj ps 2.5/ At. Now we 
can immediately evaluate the resulting linear curvature perturbation by integrating Eq. (54) to 
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get 



./bend 17 



(61) 



~ —5s(t e k-end)- 




- sin(ujAt)5s(t ek - end ) 



(62) 



(63) 



Thus, the amplitude is very similar in magnitude to the value estimated above for the process of 
ekpyrotic conversion. 

We are now in a position to calculate the variance of the generated curvature perturbation, 
which, on account of (49), is given by 



Hence, the amplitude is in agreement with the current WMAP bounds of A^(0.002Mpc _1 ) = 
(2.4 ± 0.2) x 10 -9 [9], as long as \V e k-end\ ~ (10~ 2 Mpz) 4 , i.e. the minimum of the potential has to 
be roughly at the grand unified scale for models using kinetic conversion [26]. This scale is also 
the natural scale of Hofava-Witten theory, and thus it is the scale where one would expect the 
potential to turn around. Note that for models using ekpyrotic conversion, this result implies that 
the bending must occur at a specific time, namely when the potential reaches the grand unified 
scale. In the latter models, this may or may not also correspond to the bottom of the potential. 

Finally, we should state an important assumption that we have been making implicitly up 
to now: namely, we assumed that the curvature perturbation passes through the big crunch/big 
bang transition essentially unchanged. The reason for doing so is that the perturbations we are 
considering are vastly larger than the horizon size around the time of the crunch, and hence, 
due to causality, it seems reasonable to assume that long-wavelength modes suffer no change - this 
viewpoint is discussed in much more detail in [33]. In new ekpyrotic models, in which the bounce is 
smooth and describable entirely within a 4-dimensional effective theory, this assumption certainly 
holds true. In the case of a classically singular bounce, this remains an assumption subject to 
possible revision in the future 5 . 

5 In this context, we can also mention the possibility that no conversion of entropy to curvature perturbations 
might occur before the big crunch, but that this conversion could happen during the phase shortly following the 
bang through modulated reheating [45]: if massive matter fields are produced copiously at the brane collision and 
dominate the energy density immediately after the bang, and if, furthermore, these fields couple to ordinary matter 
via a function of Ss, then their decay into ordinary matter will occur at slightly different times depending on the 
value of 5s. In this way, the ordinary matter perturbations would also inherit the entropic perturbation spectrum. 




(64) 
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IV. HIGHER-ORDER PERTURBATIONS AND PREDICTIONS FOR 

NON-GAUSSIANITY 

A. Definitions and Local Non-Gaussianity 

Now that we have seen in detail how the ekpyrotic phase generates linear, nearly scale-invariant 
density perturbations via the entropic mechanism, we can inquire as to whether the higher-order 
corrections might lead to an observable signal. We will only calculate non-gaussian corrections for 
perturbations generated via the entropic mechanism, because, as discussed in the previous section, 
this is the only robust and well-understood mechanism to date that generates ekpyrotic pertur- 
bations in agreement with observations. As we saw earlier, the linear perturbations are related 
to observations of the 2-point correlation function. Similarly, quadratic and cubic corrections to 
these perturbations are related to observations of the 3- and 4-point functions respectively. For an 
exactly gaussian probability distribution, all n-point functions for which n is odd vanish, while for 
n even, the ro-point functions are related to the 2-point function. Thus, the simplest way in which 
we could detect a departure from exact gaussianity would be due to the presence of a non- vanishing 
3-point function. 

In momentum space, the 3-point function corresponds to a configuration of 3 momenta, which 
form a closed triangle due to momentum conservation. Hence, the 3-point function is specified 
not only by its magnitude on different scales, but also by its magnitude for different shapes of 
the triangle. Or, turning this reasoning around, when we make predictions for non-gaussianity, we 
must predict both the amplitude and the shape of the momentum space triangle that we would 
like to observe. Let us make all of this more precise now. Earlier, we defined the power spectrum 
as the Fourier transform of the 2-point function, 

(C kl Ck 2 ) = (2vr) 3 ,5 3 (k 1 + k 2 )P(fci). (65) 

Similarly, the bispectrum, which is the Fourier transform of the 3-point function, is given by 

(C kl Ck 2 Ck 3 ) = (2vr) 3 <5 3 (k 1 + k 2 + k 3 )B(h, k 2 , k 3 ), (66) 

the trispectrum, the Fourier transform of the 4-point function, via 

(Ck 1 Ck 2 Ck 3 Ck 4 ) c = (27r) 3 5 3 (k! + k 2 + k 3 + k 4 )T(fci, k 2 , k 3 , fe 4 ), (67) 

and so on. The ^-functions result from momentum conservation, while B and T are shape functions 
(for a triangle and a quadrangle respectively). In the last expression, the subscript c indicates that 
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we only need to consider the connected part of the 4-point function, i.e. the part that is not 
captured by products of 2-point functions. 

For non-gaussianity of the so-called local form, it is useful to define (in real space) the following 
expansion of the curvature perturbation on uniform energy density surfaces 

C = Cl + I/nlCI + ^9nlCI (68) 

with Ql being the linear, gaussian part of £. The factors of 3/5 are a historical accident; they arose 
because this type of expansion was first defined for a different variable. In momentum space, B is 
then given by 

B = -fNL[P(ki)P(k 2 ) + 2 permutations], (69) 
5 

as can be verified straightforwardly by combining Eqs. (65), (66) and (68). Similarly, the momen- 
tum space 4-point function corresponding to non-gaussianity of the local form can be expressed 
as 

54 

T = T NL [P(k n )P(k 3 )P(h) + 11 perms.] + — g NL [P(k 2 )P(h)P(k 4 ) + 3 perms.], (70) 

where tnl and qnl parameterize the two relevant shape functions, see for example [46] for more 
details. For cosmological models in which the perturbations originate from the fluctuations of a 
single field (in our case the entropy field), tnl is directly related to the square of /nl, explicitly 

T NL = T^fNL- ( 71 ) 

Concentrating now on the bispectrum, we can see that, since for a scale-invariant spectrum P(k) ~ 
k~ 3 , we have 

1 2 2 3 3 1 

= fNL^. (73) 

This is the typical momentum dependence for local non-gaussianity [47], which is also the relevant 
one for ekpyrotic models, as we will show shortly. The signal is largest when one of the momenta 
is very small - this automatically requires the other two momenta to be almost equal, and hence 
the local form of non-gaussianity corresponds to having the largest signal generated for squeezed 
triangles in momentum space. 
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It is instructive to calculate explicitly the tree-level 3-point function for the entropy perturbation 
generated during the ekpyrotic phase, i.e. before the conversion to curvature perturbations 6 . 
Maldacena described in [48] how the expectation value for the 3-point function is given by 

((^) 3 > = -i f dt([(8s) 3 (t),H int (l/)]), (74) 

where Hi nt {t') = V sss (5s) 3 /3! = — y/en^/ {3H" 2 ) is the cubic interaction Hamiltonian. In Fourier 
space, this can be rewritten as [12, 49] 

{{5sf) = (2^) 3 <5(S i k i ) 

x2Re{-i<5s kl (t)<5a k2 (*) < 5s ka (t) / ^^ ^(O^iOK^)}- ( 75 ) 

J—oo t 

For e large, the mode functions are given approximately by {cf. Eq. (28)) 

isM ~ Wk e " U(1 ~ H»- (76) 



so that we get 

/(X a \3\ — (o^xc^.^.V. 

et 4 mr 



{{5s?) = {2nf5(^k t )^^i, (77) 



where we have used the result that 



i{\ + iht){l + ik 2 t){\ + ik 3 t)e~ iKt t J (1 - " ^ 2 t')(l - ik 3 t')e iKt ' {t')~ 5 + c.c. 

= l^ikf + ---, (78) 

with K = k\+k2+k 3 and where the dots indicate terms suppressed by powers of kit. The calculation 
shows two things: first, the momentum dependence in (77) is of the local form, and, comparing 
with (69) and using P& s {k) = l/(2i 2 fc 3 ) from (31) and (34), it corresponds to having an entropy 
perturbation of the form 

5s = 5s L + ^5sl (79) 

where 5s l is the linear, gaussian part of the entropy perturbation 5s. Secondly, the end-result is 
dominated by the terms for which \kit\ <C 1 - in other words, the dominant contribution to non- 
gaussianity is generated by the (classical) evolution on super-horizon scales, and the same holds true 
for the 4-point function also. Hence, in determining the predictions for the non-gaussian curvature 



6 Note that quantum corrections from loop diagrams will be suppressed by factors of h. 
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perturbation, we can simply perform the calculation using the higher-order classical equations of 
motion on large scales, up to the required order in perturbation theory. 

The strategy for calculating the non-linearity parameters defined in (68) is straightforward: first 
we solve the equation of motion for the entropy perturbation up to third order in perturbation 
theory. This allows us to integrate the equation of motion for £, at the first three orders in 
perturbation theory, and then we obtain the non-linearity parameters by evaluating 

rtconv-end a(2)' rtconv-end /-(3)' 

, = 5 h ek _ beg s = 25 Jt ek _ beg s 

3 / rtconv-end f(l)'\2 ' L 9 / f *cor«i-en<i /-(l)' \3 ' ^ ' 

^Jtek-beg ' ' ^tek-beg ' ' 

where the integrals are evaluated from the time that the ekpyrotic phase begins until the conversion 
phase has ended and £ has evolved to a constant value. A note on notation: we expand the entropy 
perturbation (and similarly the curvature perturbation) as 

Ss = 6sW + 6sW + Ss^ (81) 

without factorial factors and where we use Ss^ and 5s l interchangeably (in section III, where we 
dealt exclusively with linear perturbations, we sometimes wrote 5s instead of Ss^\ but hopefully 
this will not confuse the reader). 

The relevant equations of third order cosmological perturbation theory with multiple scalar 
fields was developed in [50], and we will use the results of that paper. The derivations of the 
equations are lengthy, and do not provide any further insight into the topic of this review. Hence, 
we will simply use the equations as we need them. Interestingly, it turns out that all results 
regarding the conversion process can be well approximated by various simple techniques that we 
will present further below, both for ekpyrotic and for kinetic conversion. 

B. Generation of Entropy Perturbations 

During the ekpyrotic phase, the equation of motion for the entropy perturbation, up to third 
order in perturbation theory, is given by [50] 

Ss + SHSs + V ss 5s + ^V sss ($s) 2 



+%(^) 3 + ( ¥ + 3^ + 3^ - 2^ ) (5sf5s 



+ (~^V ssa - - SH^j 5s(Ss) 2 + (lv ssss + 2^| j (5s) 3 = . (82) 



27 



Using the following expressions, valid during the ekpyrotic phase, 

V2 1 

<J = -=r V = ^ 

y/et et 2 

V2 2 2\/2e 

v -o v — — — y y 

it is not difficult to find by iteration that the entropy perturbation, to leading order in 1/e, is given 
by 

6s = 6sl + ^Ssl + e(S + J - -)<5 S |, (83) 
8 L v 60 80 5 7 L v ; 

where, just as before, 5s l oc 1/t. Note that the quadratic term agrees exactly with (79). At each 
order in perturbation theory, the non-linear corrections depend on the parameters of the potential 
at that order, cf. the parametrization of the potential in Eq. (21). The above equation specifies 
the initial conditions for the start of the conversion phase. 

Both during the ekpyrotic phase and during the conversion process, when the background 
trajectory bends, the entropy perturbation sources the curvature perturbation on large scales 
according to [50] 

2 ff 1 12 

C « ^[V, s 5s--V a 5s5s + {-V ss + - 2 V*){5 S ) 2 ] 

+ W [( "6^ " h V >°° ~ ys V M(SsfSs + (W^ + ^ 2 V S V SS + ^ s )(5 S n (84) 

The inelegant, but sure-fire thing to do now is to simply integrate this equation numerically during 
the two phases of generation and conversion, and to deduce the non-linearity parameters using 
(80), as was done in [50]. However, this approach does not give much insight into the final result. 
This is why we will present more physically transparent techniques first, which allow us to estimate 
the non-linearity parameters /jvl and g^L pretty well, and subsequently we will compare these 
estimates with the results of numerical integration and discuss the predictions. 



C. Ekpyrotic Conversion 

We start by analyzing the case where the conversion of entropy into curvature modes occurs 
during the ekpyrotic phase. For this scenario, it was shown by Koyama et al. [49] that the so- 
called 5N formalism [51, 52] is well suited. For ekpyrotic conversion, the calculation is most easily 
performed, and the result most easily expressed, in terms of the potential (18), which is why we are 
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adopting this restricted form here. In working with a parameterized potential like (21), the bending 
of the trajectory can be more complicated, in the sense that there can be multiple turns, and one 
has to decide when to stop the evolution. In this case, the results are strongly cut-off dependent, 
and without a precisely defined model specifying the subsequent evolution, it is impossible to make 
any generic predictions. 

The 5N formalism is particularly well suited to the case of ekpyrotic conversion, because the 
background evolution is simple. In fact, it turns out that by making the approximation that 
the bending is instantaneous, it is very easy to find an approximate formula for the non-linearity 
parameters at any chosen order in perturbation theory. The following calculation was first presented 
in [49] for the bispectrum, and trivially extended to the trispectrum in [50]. 

In order to implement the 5N formalism, we have to calculate the integrated expansion N = 
f Hdt along the background trajectory. Initially, the trajectory is close to the scaling solution (22). 
Then, we assume that at a fixed field value 5sb away from the ridge, the trajectory instantly bends 
and rolls off along the direction. At this point, the trajectory follows the single-field evolution 

2 2 

a(t) = (-t) 2 / c 2 02 = — ln(-t) + constant 4>i = constant. (85) 

C2 

Approximating the bending as instantaneous, it is easy to evaluate the integrated expansion, with 
the result that 

2 

N = — k In \Hb \ + constant, (86) 
c i 

where Hb denotes the Hubble parameter at the instant that the bending occurs. Note that all 
C2-dependence has canceled out of the formula above. At the end of the conversion process, we 
are interested in evaluating the curvature perturbation on a surface of constant energy density. 
But, on a surface of constant energy density, the curvature perturbation is equal to a perturbation 
in the integrated expansion [53]. If we assume that the integrated expansion depends on a single 
variable a, we can write 

( = 5N = N a Sa + ^N aa (5af + ^N aaa (5af. (87) 

In our example, we indeed expect a change in N to depend solely on a change in the initial value 
of the entropy perturbation 5s. Now, from Eq. (83), we know that 5s l oc 1/t oc H, and hence we 
can parameterize different initial values of the entropy perturbation by writing 



5s l = aH. 



(88) 
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Note that since 5s l is gaussian, so is a. With this identification, we have 

Ss = aH+ + e(g + | - j^) 3 , (89) 



so that at the fixed value 5s = 5s b, we have 



Now we can immediately evaluate 



and similarly 



a oc (90) 



qa 2 qcr 



In this way, with very little work, we can estimate the non-linearity parameters 

TNL = |/J L = \4 (94) 

25./V QaQ 25 4 . . 

ffJVL = ^iivT = l^ Cl - (95) 

We are now in a position to compare these estimates to the numerical results obtained in [50] by 
solving and integrating the equations of motion (82) and (84) 7 . There, the initial conditions were 
chosen such that they are given by the scaling solution (22), except for an 0.1 percent increase in 
the initial field velocity \^\. This causes the trajectory to eventually roll off in the fa direction, 
and to quickly approach the single-field solution (85). The results for several values of c\ and C2 
are shown table I, alongside the values estimated by the <5iV formulae. 

It is immediately apparent that the general trend is accurately captured by the 5N formulae. 
However, one may notice that the agreement is slightly less good at third order than at second, 
and also, that the 5N formulae tend to slightly over-estimate tnl an d slightly under-estimate qnl- 
But given the quickness of the 5N calculation and the complexity of the third order equations of 
motion, the agreement is pretty impressive. Of course, the 5N formula was derived subject to the 
instantaneous bending approximation. Without this approximation, we would expect a numerical 
scheme that uses the 5N formalism to yield results in close agreement with the numerical results. 



7 Approximate analytic solutions to the equations of motion were first presented in [54] . 
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Cl 


C2 


Jnl.sn 


1~NL,SN 


9NL.8N 


Inl 


TNL 


9nl 


10 


10 


-41.67 


2500 


2315 


-39.95 


2298 


2591 


10 


15 


-41.67 


2500 


2315 


-40.45 


2356 


2813 


10 


20 


-41.67 


2500 


2315 


-40.62 


2377 


3030 


15 


10 


-93.75 


12660 


11720 


-91.01 


11930 


13100 


15 


15 


-93.75 


12660 


11720 


-92.11 


12220 


13830 


15 


20 


-93.75 


12660 


11720 


-92.49 


12320 


14440 


20 


10 


-166.7 


40000 


37040 


-162.5 


38020 


41320 


20 


15 


-166.7 


40000 


37040 


-164.4 


38930 


43170 


20 


20 


-166.7 


40000 


37040 


-165.1 


39240 


44490 



TABLE I: Ekpyrotic conversion: the values of the non-linearity parameters estimated by the 5N formalism 
(using the instantaneous bending approximation) compared to the numerical results obtained by directly 
integrating the equations of motion, for potentials of the form V = —V\e~ Cl 't' 1 — V2e _C2 * 2 . 

D. Kinetic Conversion 



As discussed in more detail at the end of section III, in models motivated by M-theory a bend 
in the trajectory happens naturally just before the big crunch, during the final approach of the two 
end-of-the-world branes. This bend takes place after the ekpyrotic phase has come to an end and 
while the evolution of the universe is dominated by the kinetic energies of the scalar fields - see also 
Fig. 5. Again, there exists a simple approach to estimate the non-gaussianity parameters of the 
curvature perturbation generated by this process of kinetic conversion, and we will review it here. 
This simplified approach was first presented in [55], based on previous work in [44, 50, 54, 56]. 

This estimating procedure is based on the fact that the physics of the kinetic phase is really quite 
simple, and moreover, except for the fact that its initial conditions involve the entropy perturbation 
5s, the kinetic phase has no memory of the details of the ekpyrotic phase. In particular, only the 
total 5s in (83) matters, and the way we choose to decompose it into linear, second- and third-order 
parts is irrelevant at this point. What is more, since 5sl <C 1, the second and third order terms 
in (83) are highly subdominant and we can approximate the evolution of 5s by that of the linear 
term 5s l throughout the kinetic phase. This realization is the first ingredient of the calculation. 

The second is a compact and very useful expression for the evolution of the curvature pertur- 
bation £ on large scales and on surfaces of constant energy density [53, 54]: 

where 5V = V(t, x l ) — V(t) and a bar denotes a background quantity. This equation is exact in the 
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limit where spatial gradients can be neglected, and can thus be expanded up to the desired order 
in perturbation theory if required. If expanded, it corresponds precisely to Eq. (84) [55]. 

The third and last ingredient of the calculation is the simple relationship between 5V and 
5s during the conversion process. As we saw earlier, during the ekpyrotic phase, the curvature 
perturbation picks up a blue spectrum and is hence completely negligible on large scales. To be 
precise, since 5V ^ during ekpyrosis, there is already some conversion of entropy into curvature 
perturbations occurring at this stage. However, this contribution is entirely negligible compared 
to the subsequent conversion (note that since V s = during ekpyrosis, 5V starts out at subleading 
order), and hence we can take ((tek-end) ~ 0. Moreover, as we saw when we were discussing the 
linear perturbations, at the end of the conversion process ( is still significantly smaller than 5s, 
being given by 

Cl ~ \ss L . (97) 
5 

Hence, during the conversion process, we can take the potential to depend only on 5s. And since 
the repulsive potential is monotonic and we are interested in small departures 5s <C 1 from the 
background trajectory, it is intuitively clear that 5V is directly proportional to 5s during the 
bending. A numerical calculation readily confirms this simple relationship. 8 

During the conversion, the effect of the repulsive potential is to cause the entropy perturbation 
to behave approximately sinusoidally, as shown in Eq. (60). As we will confirm below, during the 

conversion process 5V <C a 2 , so that Eq. (96) simplifies further to 

2 Pf 

C « -pSV, (98) 

which, when integrated and upon use of (58) immediately reproduces Eq. (97). But, as argued 
above, 5s as a whole must behave approximately in this way during the conversion phase, and 
subsequently analogous relationships must hold at higher orders too: 

(2) 1 Kty/e 2 (3) 1 ,«4 4 _ \ §s 3 m 

s 5 8 L ' S 5 V 60 80 5> L K J 

These expressions immediately allow us to calculate the non-linearity parameters 

5Cf) 
3 Cl 



Inl = ^— ~ 3/t 3 Ve, (100) 

L 

9 Cl '"'60 ' 80 5' 



wlS fg!U T0(S+ !_!), (101) 



For completeness, we mention that the curvature perturbation can also be sourced by perturbations in the comoving 
energy density [41] . However, the ekpyrotic phase massively suppresses comoving energy density perturbations on 
large scales; since the kinetic phase is relatively short, they do not have time to grow and become significant, so 
that we can neglect them in our calculation - see also [56] . 
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Thus, without much work at all, we find these simple estimates for the non-linearity parameters. 

Before discussing this result, let us briefly pause to verify the approximation made in obtaining 
Eq. (98): during the kinetic phase, we can rewrite (96) as 

* = ™- < 102 » 

The approximation made above consists in writing ( pa t5V and this leads to £ pa ^tf end SV (tbend) ■ 
But we know that by the end of the conversion process £ « |<5s and hence we find that 

3t 2 bend 5V « «5s « 1, (103) 

which shows that the approximation is self-consistent and confirms the validity of (98). 

The results above indicate that the non-linearity that was present in the entropy perturba- 
tion gets transferred straightforwardly (i.e. only with an overall numerical coefficient) to the non- 
linearity in the curvature perturbation, essentially due to the simplicity of the kinetic phase. This 
calculation therefore leads us to expect no significant additional constant terms in /tvl or constants 
and ^-dependent terms in gNL', we will see shortly that this expectation is indeed borne out. 

We can now compare these estimates to the results of brute- force numerical integration [50]. 
As discussed in section III, the repulsive potential can in principle be calculated, given a specific 
matter configuration on the negative-tension brane [43]. Here, we consider four different examples 
for the repulsive potential 

V rep (x (f>2 2 , fa 2 + (P2 6 , (sinh0 2 )- 2 , (sinh^r 2 + (sinh0 2 ) -4 , (104) 

with the overall magnitude adjusted in order to obtain various values for the duration of the 
conversion (see below). These potential forms should give an indication of the range of values that 
one can expect the non-linearity parameters to take. Note that, without loss of generality, we take 
the boundary to be located at 4>2 = 0, and we only consider conversions during which 6 > - other 
cases can be related to these by an appropriate change of coordinates. 

An important parameter turns out to be the duration of the conversion [50, 56], measured by the 
number of Hubble times during which most, say 90 percent, of the conversion takes place, i.e. the 
duration is measured by the number of e-folds by which the scale factor shrinks during conversion. 
Conversions lasting less than about 0.2 Hubble times correspond to what we call sharp conversions, 
while those that last on the order of 1 Hubble time correspond to smooth conversions. It turns out 
that the estimating procedure presented above works best for the case of smooth conversions. For 
/jvl, the range of predicted values is considerably narrower as the conversion becomes smoother 
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FIG. 6: /jvi as a function of the duration of conversion and for the values k 3 — —5,0,5 and e = 36. In 
each case, we have plotted the results for four different reflection potentials, with the simplest potentials 
{4>2 2 > (sinh 4>2)~ 2 ) indicated by solid lines, while the dashed ((sinh(/>2)~ 2 +(sinh(/)2) -4 ) and dotted (r/>2~ 2 +(/>2~ 6 ) 
lines give an indication of the range of values that can be expected. As the conversions become smoother, the 
predicted range of values narrows, and smooth conversions lead to a natural range of about —50 < /jvl ^$ +60 
or so. 



[44, 56], as illustrated in Fig. 6. Moreover, /tvl can be well fitted by 

/jvx = ^3\A+5, (105) 

in good agreement with the simple estimate (100). 

We will show the results for g^i slightly more explicitly; they are illustrated in Figs. 7 and 8. In 
each case, we have plotted the results obtained for the four repulsive potentials (104). The left panel 
of Fig. 7 shows that, even more so than for /nl, the range of predicted values for g^i narrows 
drastically as the conversion process becomes smoother. In fact, for sharp conversions, typical 
values are very large in magnitude, and we expect these to be observationally disfavored shortly, if 
they are not ruled out already. Thus, phenomenologically speaking, it is much more interesting to 
focus on smooth conversions. The right panel in Fig. 7 shows that g^L is proportional to e, while 
Fig. 8 indicates that g^L depends linearly on and approximately quadratically on K3. In fact, 
all numerical results indicate that g^L scales with e, K3, K4 exactly as the third-order coefficient in 
the expression (83) for the entropy perturbation during the ekpyrotic phase, and the data can be 
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FIG. 7: On the left, gNL is shown as a function of the duration of conversion, with K3 = K4 = and for 
four different reflection potentials, and with the same line style assignments as in Fig. 6. As the conversions 
become smoother, the predicted range of values narrows considerably, allowing rather definite predictions. 
On the right, gNL can be seen to be proportional to e, i.e. proportional to the equation of state w e k- 




FIG. 8: The left figure indicates that gNL depends linearly in K4, the parameter we are using to specify the 
fourth derivative of the ekpyrotic potential with respect to the entropic direction. Similarly, in the right 
figure gNL can be seen to depend approximately quadratically on K3, the parameter determining the third 
derivative of the ekpyrotic potential with respect to the entropic direction. 



fitted by the approximate formula 



g NL ps 100 (— + — - -)e, (106) 
y v 60 80 5 J K ' 



again in good agreement with the estimate (101). 
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V. DISCUSSION OF THE RESULTS 

In this review, we have focussed on the predictions for non-gaussianity generated via the entropic 
mechanism in ekpyrotic models and their cyclic extensions. The reason for concentrating on the 
entropic mechanism is that this is currently the most robust, best-motivated and best-understood 
mechanism for generating nearly scale- invariant curvature perturbations during a contracting phase 
of the universe. As indicated by its name, the entropic mechanism achieves this by first generat- 
ing nearly scale-invariant entropy perturbations, which are subsequently converted into curvature 
perturbations. There are essentially two distinct possibilities for when this conversion can happen: 
either directly during the ekpyrotic phase (ekpyrotic conversion) or during the subsequent phase 
of scalar field kinetic energy domination right before the big crunch/big bang transition (kinetic 
conversion). In both cases, the curvature perturbations pick up non-gaussian corrections of the 
local form, although the magnitude of the second and third order non-linearity parameters Jnl 
and qnl differ substantially for the two modes of conversion. For convenience, we will repeat the 
predictions here: 9 

5 

/nl = — — c i Ekpyrotic Conversion (107) 

25 

9NL = m 4, (108) 
3 

Jnl = 2 K3v ^+5 Kinetic Conversion (109) 

9NL = (|k 4 +^«§-40)€. (110) 

For ekpyrotic conversion, the results are presented in terms of the parameters of the ekpyrotic 
potential as expressed in (18). The first thing to note is that the signs are unambiguously fixed: 
Inl is predicted to be always negative, while g^L is always positive. In order for the power 
spectrum of the perturbations to be in agreement with observations, one needs c\ > 10, which 
implies Jnl ^ — 40 and g^L ^ 2500. The current observational bounds are that Jnl = 38 ± 21, 
where the errors are quoted at la [57], while currently no strong constraints exist yet on gNL- 
These values put the predicted values for f^L for ekpyrotic conversion at 4cr or more from the 



We should caution the reader that what we are presenting here are the predictions from the ekpyrotic phase alone, 
i.e. we have assumed that the primordial density perturbations generated by the ekpyrotic phase have not been 
modified by the dynamics of the big crunch/big bang transition, nor that there are additional effects of relevance 
during the first stages of the subsequent expanding phase of the universe (it is conceivable, for example, that 
models might be constructed which utilize the ekpyrotic phase, but where there are additional fields that are 
relevant during the expanding phase and which would contribute to the curvature perturbation - in this case, the 
final predictions will evidently be model-dependent). 
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central value, and thus this type of conversion is on the verge of being ruled out by observations. 
What this means for model-building is that models that have made use of the ekpyrotic conversion 
mechanism, such as new ekpyrotic models [10, 12], might have to be modified in a way such that 
the conversion occurs only once the equation of state becomes small. This is not unnatural, as the 
ekpyrotic phase must come to an end anyhow before the bounce phase begins, and if the conversion 
occurs after the end of the ekpyrotic phase, the predictions will be closer to those predicted by the 
kinetic conversion mechanism, which we turn our attention to now. 

For kinetic conversion, the results (109) and (110) are presented in terms of the parameters of 
the ekpyrotic potential as written in (21). The parameters k^,k,4 are expected to be of 0(1). The 
fast-roll parameter e is typically of O(10 2 ), with a lower bound of about 50 in order for the power 
spectrum to be in agreement with observations. Thus, the fitting formulae predict the bispectrum 
parameter /nl to be of order a few tens, with the sign being typically determined by the sign 
of K3. This is in good agreement with current bounds for very natural values of the parameters, 
such as e = 100 and — 1 < K3 < 5, for example. The associated prediction for the trispectrum 
qnl is rather interesting, owing to the fact that the K3, ^-independent contribution is quite large. 
It sets the "central" value of the prediction at — 40 e and thus tends to make qnl of order a few 
thousand and negative in sign. Even with largish values of K3 ~ 5 and K4 ~ 5, the prediction is 
still negative g^L ~ — e> and hence the negative sign of g^L is a rather robust prediction of the 
kinetic conversion mechanism. For completeness, we note that the second trispectrum parameter 
tml is always given in terms of Jnl according to (71), and thus it is predicted to be positive and 
of order a few hundred. Observational limits on tnl will thus provide a consistency check of the 
predictions, and will give an indication whether or not the fluctuations of a single field (in this case 
the entropy field) were at the origin of the final spectrum of curvature perturbations. 

What are the implications of these predictions? Assuming that the observational error bars 
quoted above will shrink in the near future and that the ekpyrotic conversion mechanism will 
be ruled out, we will focus here on the predictions of the kinetic conversion mechanism. The 
most interesting feature is that the natural ranges of the non-linearity parameters are at a level 
that will be observable by near-future experiments, perhaps already with the PLANCK satellite. 
Hence, in a few years, we will know whether or not ekpyrotic models, in combination with the 
entropic mechanism, will be viable or even preferred by the data. In this vein, it is useful to briefly 
contrast the predictions discussed here with those of inflationary models 10 : for simple, single-field 



Here, we only contrast the predictions for non-gaussianity of the local form. Non-canonical kinetic terms addition- 
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inflationary models the predicted values for all non-linearity parameters are very small, of order 
1 or smaller [48, 59]. Even though these values lie within the predicted ranges of the kinetic 
conversion mechanism, it is clear that, since these values arise naturally for single-field inflationary 
models whereas obtaining the same small values in ekpyrosis would require a very finely tuned 
potential, in case of a non-detection of non-gaussianity the simple inflationary models will be 
strongly preferred. Conversely, since multi-field inflationary models can give rise to just about any 
values of the non-linearity parameters [60], in case of a detection of non-gaussianity compatible 
with the values predicted by the kinetic conversion mechanism, those latter models will be preferred 
by the data. At that point it will be necessary to include also the observational evidence for or 
against primordial scale- invariant gravitational waves. Such gravity waves would strongly favor 
inflationary models, whereas their absence is compatible with ekpyrotic models. The absence of 
primordial scale- invariant gravitational waves, combined with local non-gaussianity parameters in 
the ranges predicted by (109) and (110), would be nothing short of spectacular, as they would point 
towards the existence of a contracting phase prior to the big bang, and open up the possibility 
of a multiverse in time, with correspondingly vast timescales unlike anything we are used to in 
cosmology right now. 
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